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O ■ Abstract 

<N ■ 

f^ ' The longitudinal susceptibility xl of the 0{N) theory in the broken phase is ana- 

^H' lyzed by means of three different approaches, namely the leading contribution of the 

^^ ' 1/A^ expansion, the Functional Renormalization Group flow in the Local Potential 

approximation and the improved effective potential via the Callan-Symanzik 
equations, properly extended to d = 4 dimensions through the expansion in 
powers of € = 4: — d. The findings of the three approaches are compared and their 
agreement in the large A^ limit is shown. The numerical analysis of the Functional 



^' -1 

D ■ Renormalization Group flow equations at small N supports the vanishing of Xl in 

d = 3 and d = 3.5 but is not conclusive in d = 4, where we have to resort to the 

Callan-Smanzik approach. At finite N as well as in the limit N — t- do, we find that 

^ • Xl vanishes with J as J*^'^ for e > and as (ln( J))~^ in d = 4. 

m . I. INTRODUCTION 

^' 

^^ . The 0{N) scalar theory is the simplest field theory with a continuous global symmetry, 

hence it has been widely used as a test to check new methods and analytic or numerical 
approximation schemes and its properties have been studied in detail. A very important 

^ I property that is realized in this theory is the spontaneous symmetry breaking of the global 
symmetry, a phenomenon that plays a central role both in quantum field theory and in the 
study of phase transitions in statistical physics. In particular, the 0{N) theory is the simplest 
model suitable for studying the universal properties of the QCD chiral phase transition at 
finite temperature, as well as other aspects of the chirally broken phase of quark models. 

The spontaneous breaking of a global symmetry is signaled by a nonzero value of the order 
parameter, which in the 0{N) theory corresponds to a nonvanishing vacuum expectation value 
of one component of the scalar field, and is typically associated with the appearance of massless 
excitations, the Goldstone bosons[l|, whose propagation induce infrared divergences that can 
substantially modify the mean field indications. The nature and the strength of the infrared 
divergences is strictly related to the number of Euclidean dimensions d of the system and in 
particular the Coleman-Mermin- Wagner theorem, specifically formulated for spin models in|2| 



and for quantum fields in |3| , precludes the existence of a phase with conventional long range 
order, or spontaneous symmetry breaking in c? < 2. When d increases, the infrared divergences 
become less severe and one has to resort to a suitable approximation scheme in order to deal 
with the infrared sector and extract finite answers for the physical observables. 

A specific observable of the 0{N) theory where the effect of the Goldstone bosons is crucial, 
is the mass of the field that acquires a nonvanishing vacuum expectation value ipo, indicated as 
the longitudinal field. In fact, while it is possible to make use of general symmetry arguments to 
put constraints on the mass of the Goldstone bosons associated to the other (A^ — 1) transverse 
fields, the same is not true for the longitudinal field. When the scalar theory is coupled to a 
magnetic field J, the Ward-Takahashi identities, which relate the various correlation functions 
of the theory according to the symmetry properties of the problem, state that the inverse 
transverse susceptibility, defined as the inverse transverse propagator at zero momentum, is 
Xx^ = ^t{p = 0) = J/(po- Then, when the magnetic field J is turned off the Ward-Takahashi 
identities ensure the vanishing of Tt{p = 0), that is a vanishing Goldstone bosons mass j^. 

There is no such constraint on the inverse longitudinal susceptibility Xl^ = ^l{p = 0) and 
therefore on the mass of the longitudinal field. However, the development of Renormalization 
Group techniques provided a deeper understanding of the infrared structure of the theory, 
and for this specific problem in 2 < (i < 4, the relation Xl^ °^ j(4-'^)/2 jg predicted (5|-[7|. 
This scaling of the longitudinal mass has been repeatedly analyzed and confirmed in various 
approaches |8l-ll2l|. However the behavior of this quantity in d = 4 is different because in 



this case the power law scaling is replaced by much weaker logarithmic corrections. In 
a few arguments are presented to show that xZ^ ^-Iso vanishes in d = 4 and in [l3|, where 
the two point function of the 0{N) theory is studied in the framework of the Functional 
Renormalization Group, a numerical investigation suggests a vanishing mass in c? = 4. 

The Functional Renormalization Group is a reformulation of the Wilsonian Renormalization 



Group [IJ], based on the infinitesimal integration of momentum modes from a path integral 
representation of the theory with the help of a Wilsonian momentum cutoff. The resulting 
functional flow equations interpolate between the microscopic theory at short distances and 
the full quantum effective theory at large distances. Various realizations of the Functional 
Renormalization Group were developed and, among the most renowned, there are Polchinski's 



and Wetterich's formulations 15|, llGJ (for reviews see [17|] [l8|). Although this method is 



particularly suitable to analyze the critical domain and is especially powerful in the evaluation 
of the related critical exponents (see for instance (l9l-l2l| for an extended list of references), 
it also turns out to be an important tool in the study of nonperturbative features of the 
effective potential and effective action. For instance, the effective potential of the broken 



phase of scalar theories with its convex form (this convexity property, 22|-|24l| . does not show 



up in perturbation theory 25|] and can only be recovered in some nonperturbative scheme 26| ) 



has been largely investigated by means of Renormalization Group techniques 27H38|], 13|. In 



particular, an improvement on the determination of the effective potential in the ordered phase 



of the 0{N) theory is considered in [13|, where a numerical computation of the longitudinal 
propagator which includes the wave function renormalization, i.e. momentum dependent 
corrections to the Local Potential approximation, is performed in d = A with the result that 
large wave function renormalization factors appear at least away from the critical region of 
the phase transition to the disordered phase. 

In this paper we reconsider the study of the longitudinal mass in the broken phase of the 
0{N) theory, in the framework of the Functional Renormalization Group, generalizing the case 



d = 4: examined in [13[ to the range 2 < c? < 4, however not including momentum dependent 
corrections that are not crucial for our purpose and therefore limiting ourselves to the Local 
Potential approximation. We first analyze the problem in the simpler large N limit case and 
compare the results with those obtained at lowest order in the 1/A^ expansion. Then we extend 
our investigation to small values of A^, also confronting the results with the Renormalization 
Group improved effective potential as obtained from the Callan-Symanzik equations. With 
this compared study, we gain analytic control on the problem which allows to overcome the 
limits of the numerical analysis, that are particularly evident in proximity of rf = 4, and finally 
we get a more complete picture of this specific issue. In Sec. [TTl we briefly go through the 
determination of the longitudinal susceptibility at the leading order in the 1/A^ expansion and 
in Sec. IIIII we discuss the Functional Renormalization Group approach in the same limit. In 
Sec. IIVI the problem at small N is investigated by means of the Functional Renormalization 
Group techniques and of the perturbative Renormalization Group improvement. Conclusions 
are reported in Sec. |Vl 

II. THE LARGE A^ LIMIT 

We start by considering the 0{N) theory defined by the partition function 

Z[J]= /"p$e"^[*l+^^-* (1) 

where $ is a A^-component scalar field and the action S[^] reads 

5[$] = |d^x{i[9,$f + ^$^ + A[$2]^|, (2) 

and we focus on the limit of large number of field components, A^ — ?■ oo. By using the standard 
technique of introducing an auxiliary field p in the partition function via the following constant 
multiplicative factor (where the integration contour in the complex p plane must be taken 
parallel to the imaginary axis) 



and then, by indicating with ip the longitudinal field component parallel to the source J and 
integrating over the remaining (A^ — 1) components of $, the new effective action, functional 
of the fields ip and p, that appears in the partition function 

Z[J]= /p^Ppe-^^-f-f'^'^l-^'^'^^^'^) (4) 



is: 



SefAv^p] = / d'^X 



1,^ ,2 1 2 3M2 3 2 



+ !:L^tT\n[{-d'+p)6''ix-y)] (5) 

where we also included the coupling to the source J. 

This procedure shows the explicit dependence on A^ and it is understood that, for large N, 
the square field v?^, the inverse coupling 1/A and the effective action Sejf are proportional to 
N. The limit A^ — )• cxo is typically studied by keeping the product (A^ A) finite and performing 
the functional integral by means of the steepest descent method. In particular, by restricting 
to uniform configurations for ip and p, one gets the two coupled extremum equations: 



6S, 



eff 



6ip 



pip = J, (6) 



- ip-M') + ^ + - j7d^-r— = (7) 



6p A'"^ ' 2 2 J {2-kYp^ + p 

and the fiuctuations around the solutions of Eqs.( l6]l71) are neglected because they are sup- 
pressed by inverse powers of A^. 

From Eqs . (1^171) at J = one can immediately find the symmetric phase solution where 
(^ = and p is the selfconsistent square mass solution of the gap equation (JTj). Clearly this is 
possible only if M^ > M^, where 

_2 AAT /• d'^p 1 

6 J (27r)^p2 V ) 

For M^ < M^ the system admits another solution and the broken phase is realized. 

The broken phase, which will be considered below, corresponds to a nonvanishing field 
(y5 7^ at J = 0. This, according to Eq. (|6]) , implies p = and at the same time Eq. ([7]) 
yields the value of the field p^ that correponds to the minimum of the potential 

2 6M2 6M2 , . 

^1 = -^^^- (9) 



Then we are able to rephrase Eqs. (1^|7|) at finite J, in terms of ipo 

(10) 



J 
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-n[ "■'" 


-1 
/ 1 
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a'' 


VP^ + P 


p2 



^'-^l = ip-N j^{-^—:-i, (11) 



and it is easy to check the dependence of the field (p on the source J, at least in 2 < c? < 4, 
where the integral in the right hand side of Eq. flTTl) is finite and proportional to p^~ . After 
performing the integral, Eq. flTTl) becomes 



,2 ,._6p ,^r(i-rf/2)^,^,^ 



^^-^o = f-iV^^^^p|p|-- (12) 

Eq. ( fT2l) can be explicitly solved for p, at least in the region \p\ << 1, (which means close 
to the minimum of the potential, according to the saddle point equations (16|7|1 with ip y^ ) 
where the term p^~ is dominant with respect to the linear term p, 

The solution (fT3l) . when replaced in Eq. (j6]), V'{ip) = ipp, yields an analytic expression of the 
derivative of the potential in the region \p\ << 1. 

Incidentally when M^ = M^, i.e. ipo = 0, one finds: 



d-2 



if oc J'i+'2 (14) 

and the inverse of the exponent of J gives the correct critical index 6 = {d + 2)/{d — 2). This 
picture is slightly modified in the broken phase with ipo ^ 0, again for 2 < rf < 4. In fact in 
this case it is convenient to define the shifted field ipsh = V^ — V^o so that in the left hand side 
of Eq. f lTT]) . for J — > one can retain the 0{(psh) term, neglecting 0{iplfj and, in the same 
way, the leading behavior of p is given by p ~ J/<Po- Then, as before, in the right hand side 
of Eq. flTTj) one can retain only the leading term proportional to p~2~ and finally one finds 

(Psh oc J^~ . (15) 

It is then straightforward to get the behavior of {dip{J)/dJ) from the derivative of Eq. (fT5l) 
for small J and 2 < d < 4, by recalling that ipo in Eq. ([9]) is independent of J, 



,4-d 



xi\j) = {^)"o,r-^. (16) 



Then, clearly, for 2 < ci < 4 when the source J is turned off we have Xl^{J = 0) = 0- 



In d = 4, the integral in Eq. fITT]) contains a logarithmic divergence and the expression 
in Eq. f lT^ is modified by logarithmic corrections so that it is no longer straightforward to 
generalize Eq. ( !T6l) to rf = 4 by following the same procedure. Therefore, we shall reconsider 
the problem by analyzing in detail the behavior of the integral, including the case d = 4. In 
particular we use Eqs.f jTOfTT]) to express p = p{(p, J) and finally y? as a function of the source 
J: ip = ip{J). Therefore, by deriving Eqs. fllOflip with respect to J and putting them together 
one has 



Xl 



\ dJ 



l + Nlt\p) 



(17) 



where 



M) 



I^^'ip) 



d'^p 



(27r)'^ (j92 + pY 



The inverse longitudinal susceptibility is given by Eq. fITTj) in the limit J — )■ 0. Incidentally 
we note that the same result as in Eq. flT7|) is obtained as well, when computing the longitudinal 
mass m\{J) = Xl^{J) from the inverse propagator of the theory SSeff/{Sip6ip) (by treating 
p as p{(p) ) at zero momentum. To extract the correct behavior of Eq. flT7|l when J — )■ 
it is essential to analyze in detail the integral /g (p) including the case d = 4. We compute 
I2 (p) by inserting an ultraviolet Euclidean cutoff A which protects the result from possible 
ultraviolet divergences: 



rid) 



irip) = c,\p 



.-4 /-^^/l^l , x^ 
dx 







Cd\p\ 2 



Cd\p\ ' 



x + lf 
{d-2) 



1 + F 



d 



P 



"'2' 2 



-^^r(i-rf/2)r(i + rf/2) + ^(M 



d 



\p\ 

4-d 



2\d/2 



6-d 

Ad /|p|^ 2 



(rf-6) \A^ 

4 (v{2-d/2) 



O 






\p\ 



(47r)^/2 ' {d-A){ATxY/'^T{d/2) \h? 



(19) 



where Cd= {{A'KYI^V{d/2))- 
Hypergeometric function. 



r(x) is the Gamma function and 2-^1 (a, &; c; z) is the Gaussian 



We note that Eq. ( 1191) when d < A does not present any ultraviolet divergences. The leading 
term is [r(2 — (i/2)/(47r)'^'^]|p|^~ and all the other terms in the development of Eq. ( iT9l) are 
suppressed by positive powers of \p\/ h^. If we take A — )■ cx), only the leading term survives 
and we end up with the same result already shown in Eq. (TT6l) . However, due to the presence 
of the cutoff A, we are able to analyze the case d = A. In fact one can easily check that, by 



6 



expanding the last line of Eq. f lT^ in powers of e = 4 — rf, in the limit (i — )■ 4 the divergent 
terms in e cancel out leaving a logarithmic term : 

In this case the ultraviolet cutoff A appears logarithmically and we have to recall that for 
d = 4 the limit of infinite A corresponds to a trivial behavior of the theory that turns out 
to be noninteracting at any finite renormalization scale. Therefore we limit ourselves to the 
study of an effective theory with fixed ultraviolet cutoff A. Then the linear term in p in Eq. 
( 1T7|) is negligible and the inverse susceptibility goes to zero as (| ln(J)|)~^ , when J — )■ 0. We 
shall analyze in more detail this aspect in Sec. IIVI 



III. FUNCTIONAL RENORMALIZATION GROUP AT LARGE N 

In this Section we shall use the Functional Renormalization Group flow equations to analyze 
the theory in the large A^ limit. According to this approach one can follow the evolution of 
the running effective action r^ as a function of a momentum k scale which, starting from 
the ultraviolet cutoff A, goes down to the infrared limit k = 0. The flow of F^ starts in the 
ultraviolet region with an initial boundary condition given by the bare action and its final 
endpoint at k = corresponds to the full effective action of the theory. The integration of the 
fluctuations down to the scale k, which generates the flow of F^ , is performed by modifying the 
action of the theory with an additional term, quadratic in the fieds, that contains a /c-dependent 
regulator Rk which acts like an infrared cutoff that effectively enables the integration of the 
modes above k. In order to fullfill the requirements that Tk=A is the bare action and Tk=o 
coincides with the full effective action, the regulator Rk must have the following asymptotic 
properties : Rk^oa = oo and Rk=Q = 0. Then if $ = ($i, . . . , $Ar), F^ (g, — g; ^)ah is the 
second functional derivative: 

' ^^' ^'^^"^.^ 5$,(g)5$,(-g)' ^''^ 

[Rk{<i)\ab is the /c-dependent regulator and t is defined as t = ln(A;/A), the flow equation reads: 

dtU^] = lj-0yi^' {^^R^i^) [rr[?> -?; ^] + Ml)] "') (22) 

where the trace is performed over the internal indices. 

In order to analyze the behavior of the longitudinal susceptibility for the 0{N) theory, it 
is sufficient to approximate Eq. f l22|) to the so called Local Potential approximation, where 
we limit the form of Ffc[$] by ignoring all possible scale dependent coefficients of the kinetic 
{0{d^)) term and by fully neglecting 0(9^) terms, 

r.[$] = ^ (Ukim + I [d,^ad,^a] + 0(8')^ . (23) 



In Eq. fl25]) . due to the symmetry of the problem, the potential Uk depends on the modulus 
|$| = V^a^a and, according to our approximation, Eq. ( 122|) reduces to a flow equation for 
the potential f/fed*^!), with ultraviolet boundary conditions fixed in Eq. ([2]) : 



f/fc=A(|$|) 






(24) 



and the final infrared limit of the flow, Uk=Q{\^\)i corresponds to (an approximation of ) the 
full effective potential. 

The explicit form of the flow equation depends on the particular choice of the matrix Rk, 
and here it is convenient to take the following regulator that allows us to solve analytically 



the momentum integrals 39 



-42 



[Rk{q)U={k^-q^) Q{e-q')6, 



ab 



(25) 



where 9(x) indicates the Heaviside function. With the regulator in Eq. (125|) . the flow equation 
for Ukim) is 



dtUkim) 



2C, 
d 



^2H-d 



(N 



+ 



1 



LP + f/Ki$i)/i$i k^ + u'^m)\ 



(26) 



where C^ is defined after Eq. fITI?]) and the primes indicate derivatives with respect to |$|. 
The contribution of the {N — 1) Goldstone bosons and that of the longitudinal field in Eq. 
( 126|) are clearly evident. 



The very simple structure of Eq. ( HE^i allows us to analyze the large N limit, where the 
contribution of the longitudinal component is negligible and therefore, by deriving with respect 

to |$|, we get 



dkU'.m) 



2C, 



<i lA+d 



N 

W\ 



um\)/m-u'>:m) 



(27) 



We recall that the inverse longitudinal susceptibility coincides with the square mass of the 
longitudinal field which is the second derivative of the effective potential evaluated at the 
vacuum expectation value of the field, i.e. at the minimum of the effective potential. Therefore 
in this framework it corresponds to ?7^'^q(|$o|) with |<l>o| such that f/(.^Q ( | $0 1 ) = . The flow 
equation in the form of Eq. (I27|) allows us to make some statements on t/^'^o ( I "^o I ) • 

Before doing that, we have first to recall the behavior of the running potential in 
the broken phase when /c ~ 0, which has repeatedly been described in various contexts 
32| . 35|-|38|, 131]. Essentially, it turns out that for k much smaller than the other 



29 , 31 , 30 



characteristic scales of the problem, the derivative of the potential is very well approximated 
by Uj^im) = -A;2|$| + 0(P, |$n in the region < |$| < |<l>o|, while for |<l>| > |$o|, f/fc(|$|) 
has already reached its infrared limit and is practically k - independent with the exception of a 
very narrow region around |$o| where the two different branches merge. In this narrow region, 
no discontinuity is observed in ?7^(|<l>|), and in particular one has f/^^Q(|$o|) = 0. Accordingly, 
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Figure 1: Left Frame. The derivative of the effective potential in the large N limit given in Eqs. 
([B]) and d?]) (dashed) compared with the Functional Renormalization Group {7^(|<l>|/\/]V)/\/]V at 
k = 0.002 (solid) for d = 3, PP = -1, NX = 2.4, A = 10. 

Right Frame. Zoom of the Left Frame in the region around ifo = 1.47046 , with the inclusion of the 
analytic solution in Eq. (fTSl) (dotted). 



one finds that f/^'(|$|) ~ k"^ for < |$| ^ |$o| ^-i^d therefore, when /c — !■ 0, its left limit 
(when |$| — |$o| — ^ 0^) vanishes. Our goal is to establish whether or not its right limit (when 
|$| — |$o| — ^ 0+) vanishes, i.e. whether or not f/^'^Q(|<l>o|) is discontinuous. 

In order to clarify this point, we observe that according to the mentioned solution for 
t/^(|$|) at k ^ 0, the denominator in Eq. (1271) . computed at |$o|) vanishes as k^ and, when 
combined with the other powers of k in Eq. ( 1271) , it gives the factor /c*^"^ which diverges in the 
infrared limit /c — )■ (in our problem d < 4). As the left hand side of Eq. (1271) at |$o| and for 
k ^ cannot be divergent, the numerator in the right hand side must vanish rapidly enough 
to compensate the diverging factor. Since the continuous function [/(.(|$o|) vanishes as k"^, it 
follows that also f/('(|$o|) niust vanish when A; — )■ and therefore no discontinuity can appear 
because its presence would produce a divergence in the right hand side of the differential 
equation when /c — )■ 0. This analysis indicates that in the large A^ limit the flow equations 
constrain the inverse susceptibility to vanish as ?7^'(|$o|) = ~^^ + 0{k^). Unfortunately this 
argument cannot be straightforwardly extended to the case of finite N, because of the presence 
of the longitudinal contribution in the differential flow equation. 

On the other hand, one can try to extract the longitudinal susceptibility by numerically 
solving the Functional Renormalization Group flow. In fact the numerical resolution of the 
flow consists in the integration of a partial differential equation for the potential ?7fc(|$|) with 
fixed ultraviolet boundary conditions indicated in Eq. (^^. The integration is carried out with 
the infrared scale that approaches k = and the corresponding output of the flow provides a 
numerical approximation of the effective potential. Further details of the numerical integration 




Figure 2: Left Frame. The inverse longitudinal susceptibility in the large A^ limit given in Eqs. (|T< 
and d?]) (dashed) compared with the Functional Renormalization Group U'^{\^\/\/N) at k = 0.002 
(solid) with the same parameters used in Fig. [H 

Right Frame. Zoom of the Left Frame in the region around ipo = 1.47046, with the inclusion of the 
analytic solution in Eq. p3]) (dotted). 



can be found in 13 . 



The numerical output of the Functional Renormalization Group flow can be compared with 
the solution of the saddle point equations obtained in the large A^ limit. In particular, we 
can compare f/(. and f/^' with the first and second derivative of the potential obtained from 
Eqs. f l6llT7|) where the parameter p is replaced according to Eq. ([7]). We can also check the 
analytic solution for p displayed in Eq. flT3|) . which is valid only in a narrow region around 
the minimum. As a specific example, in Figs. [1] and [2] the first and second derivative of the 
potential are displayed for d = 3 and with the following set of bare parameters, M^ = —1, 
A^A = 2.4. The ultraviolet cutoff is taken A = 10 for the Functional Renormalization Group 
and it is rescaled of a factor 2/3 in Eqs. fl6|71IT7|) . because of the difference in the dependence 
of the renormalized parameters on A, as discussed in 13 . 

In the limit A; — )■ 0, the Functional Renormalization Group flow produces a convex potential 
with a fiat part in the region | $ | < | $o I ) 27H38| , 13| , in agreement with a basic property of the 
effective potential 22|-|24| . Actually, the numerical routine adopted here fails to integrate the 
flow equations below a particular value of the infrared scale which corresponds to k = 2 ■ 10~^ 
in the example considered, and this explains the small difference from zero observed in the 
Functional Renormalization Group curve of Fig. [1] (right frame) for |$| < |$o|- It is also well 
known that the convexity property of the potential is not recovered when using approximations 
based on an expansion around a saddle point, such as the loop expansion or the large N 
expansion considered in Sec. [Til In these cases, it is necessary to consider the contribution of 



all the saddle points in order to obtain a convex potential 25[. This explains why in Fig. [T] 
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Figure 3: The inverse longitudinal susceptibility in the large A^ limit plotted as a function of the 
source J. All parameters are fixed as in Figs. [T]and[2j The same coding of Fig. [2] is used for dotted, 
dashed and solid curves. 

we observe that at the leading order in the 1/A^ expansion, despite its claimed exactness, the 
effective potential does not present the expected convex shape in the region between |$| = 
and |$| = |$o|. 

On the other hand, for |<l>| > |$o| the computation performed by considering a unique 
saddle point is exact in the limit N ^ oo and this explains the good agreement of the various 
curves above |$o|- We observe in the right frame of Figs. [1] and [2] that the analytic solutions 
for J (i.e. the derivative of the potential) and Xl^ (see Eqs. (PI7|17|13p ). perfectly reproduce 
the numerical solution of the large A^ equations only for |$| very close to l^ol^ as expected. 
From the plot of the curves in Fig. [2] it is also evident that the second derivative of the 
potential, i.e. the inverse susceptibility, vanishes at |$o|- Finally the plot of U'f! computed via 
the Functional Renormalization Group is in good agreement with the other curves (compatibly 
with unavoidable numerical fluctuations) despite the flow is truncated before reaching k = . 

Fig. [3] shows the inverse suceptibility as a function of the source J with J > for the same 
set of parameters of Figs. [T]and[2l Again the numerical solution of the saddle point equations 
dnj and ([7]) is displayed together with the analytic solution in Eq. ( 1T3|l and with the output 
of the flow equation. Fig. [3] shows the range of small values of J where the flow equation is 
no longer accurate. For larger values of J, the dashed and the solid lines coincide. In addition 
we notice that the analytic solution obtained by replacing p in Eq. ( IT71) with the result of 
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Figure 4: U^ plotted versus J for d = 3 and A; = 5 • 10 ^ (dotted) , d 
d = 4 and k = 2 ■ 10~^ (dashed). Values of the other parameters : M^ = 



0.05 



3.5 and k = 10^^ (solid), 
-1, A = 2.4 and A = 10. 



Eq. (I13p . and which displays the scahng behavior indicated in Eq. ( 1T6|) . does in fact coincide 
with the full numerical solution only for extremely small values of J. In this very small range, 
the Functional Renormalization Group flow equation solution has already deviated from the 
correct behavior and a more efficient algorithm becomes necessary to push k closer to zero to 
obtain a better agreement in this region and to recover the scaling behavior of Eq. ( IT6l) . 



IV. SUSCEPTIBILITY AT FINITE N 

We start the analysis of the problem at finite A^ by looking at the full fiow equation of the 
potential in Eq. f l26|) with both contributions of the transverse and longitudinal fields that 
introduce the dependence of the fiow on two different scales, namely the running masses of 
these two fields. As already noticed it is not easy to put analytic constraints on the longitudinal 
mass when fc — )■ due to the presence of these two separate contributions. 

We turn to the numerical resolution of the fiow equation to extract some physical indi- 
cations in the infrared limit. In particular, we are interested in the behavior of the inverse 
susceptibility, i.e. the mass of the longitudinal field. Let us first consider the case of A^ = 4 
and bare parameters M^ = —1, A = 2.4 and cutoff A = 10. In d = 3, the flow is stable 
down to A; = 5 ■ 10~^ and the corresponding solution U^ is displayed in Fig. H] (dotted line) 
as a function of the source J = f/p with J > 0. The other two curves in Fig. H] respectively 
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correspond to d = 3.5 (solid), where the lowest reached infrared value is /c = 10 ^, and d = 4 
(dashed), with k = 2- 10~^ and the bare parameters have the same values as in the d = 3 case. 

Due to the finite value of k, the three curves in Fig. lUdo not reach the origin (J = 0, U" = 
0). However, while for smaller d the curves clearly approach the origin, the same is not true 
when d is at the upper critical dimension. In particular, when d = 3 (dotted) the infrared cutoff 
k can be pushed very close to k = and the corresponding curve in Fig. H] smoothly bends 
toward the origin. When d = 3.5 (solid) we already observe numerical fluctuations in the plot 
but even in this case there is a change of slope and the curve approaches the origin. In d = 4 
although the infrared cutoff is not much larger than for d = 3.5, the numerical fluctuations 
are stronger and make the deep infrared region [k — )■ 0) inaccessible. The corresponding curve 
(dashed) does not show any bending toward the origin and it is not possible to establish from 
the plot whether or not the second derivative of the potential vanishes at J = and k = 0. 
In fact, as we shall see below by applying standard Renormalization Group techniques, such 
problem is related to the fact that power law corrections are replaced by smoother logarithmic 
corrections in d = 4, which turn out to be more difficult to treat numerically. For better 
understanding this point, we now proceed to the study of the theory at finite A^ and for d 
close to 4, in the framework of the Callan-Symanzik approach. 

To this end, it is convenient to rewrite the potential as 

VU) = —02 + ^ , (28) 

^^' 12^ 24A' ^ ' 

where we have introduced the rescaled field 

02 = A$2 , (29) 

and the "temperature" r 

r=jiM'-M'^), (30) 

is normalized in such a way that it vanishes at the critical value M^ = M^ that separates the 
symmetric and the broken phase. It is understood that M^ includes corrections with respect 
to Eq. ([8]) which is obtained in the A^ — > oo limit. We are specifically interested in the broken 
phase and therefore we take r < 0. 

We also recall that the coupling A has dimension A^"*^, expressed in units of the ultraviolet 
cutoff A, and the corresponding dimensionless coupling is defined as m = AA'^"^. By defining 
t = In ^ as the logarithm of the scale n, (we avoid to use the same notation of the running 
momentum k adopted for the Functional Renormalization Group fiow equations due to the 
different roles of the two scales), the /3-functions of the various parameters in the potential, 
as given by perturbation theory at lowest order, are {e = 4 — d) : 

/3 = ^ = -e^ + Cd^^u\ (31) 

d{M^-Ml) N + 2^^^^. ^. 



dt 



C,^-u{M' - M^). (32) 
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Note that $ at lowest order has no correction and therefore no t dependence. A word of 
caution concerning the usual way of obtaining these /3-functions has to be said: the fact that 
the 0{N) model in the broken phase contains more than one scale requires a more careful 
treatment, J43|, and we shall come back again to this point. From Eqs. ( 129)1301131 Il32l) we get 






-77 = -'^CdT 



dt 



e 1 du _CdN + 8 
2^2u~dt ~ ~2 



-u . 



(33) 
(34) 



By solving the differential equations related to these /3-functions, one can determine the 
evolution of the various parameters with the scale /i, after having set the boundaries at the 
bare values when /i = A. Below, the dependence of the various variables on the parameter t 
is explicitly displayed while it is omitted for the bare variables, which correspond to the value 
t = 0. So the dimensionless coupling is 



u{t) = ue 



-et 



3e ^ ' 



(35) 



while the dimensionful coupling A(t) = u{t) /i*^ reads 

1 



m 



1 GKAH 
A 3e 



(/x- - A-) 



and the temperature r is 



r(t) = r 



u 



-et 



U{t) 



6 

N+e, 



(36) 



(37) 



Finally, the renormalization factor S,{t) of the rescaled field is 

-1/2 



m 



olodt'^^it') 



U 



-et 



U{t) 



(3^ 



We are now able to evaluate the Renormalization Group improved potential Vrg/(0); so- 
lution of the Callan-Symanzik equation. This is achieved by replacing in Eq. (1251) the bare 
quantities with the corresponding scale dependent parameters: 



V, 



r(t) 



RGI\ 



m 



2i2 



-X\t)^'- 



(39) 



12 "'"' " 24A(t) 
The derivative of Vrg/(0) (in this Section the prime indicates derivative with respect to 0) is: 



vkaM) = ^em 



r{t) 



1 



m 



m 



2 j2 



J 

7f 



(40) 
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We can also determine the second derivative of Vnciit) which, after recalhng the definition of 
J in Eq. (^, reads 

Let us now go back to a somewhat dehcate point. The model under investigation has two 
characteristic scales Mi(0) and M2(0), respectively associated with the mass of the longitu- 
dinal field and with the Goldsone bosons mass: Mf (0) = XV" {(p) and M^{(f)) = XV' {(p)/ (p. If 
these two scales are comparable, then the usual procedure to compute Vrg/I^), that we have 
outlined above, is essentially correct. On the other hand, when these scales are significantly 
decoupled, new aspects, that we have overlooked up to now, appear 43|. It is important 



to observe that some perturbative diagrams contributing to the longitudinal two point func- 
tion are affected by infrared divergences due to the loops of Goldstone bosons and a reliable 
estimate of V^g7-(0) is obtained only when the Goldstone contributions are included in the 
evaluation of the Renormalization Group functions. Indeed, the Goldstone bosons have zero 
mass in absence of external source and, when J is turned on, their square mass is proportional 
to J, while the longitudinal field starts with a nonvanishing bare mass at t = and, at least 
in the A^ — )■ oo case, its renormalized square mass is proportional to J"^/^ when J approaches 
zero. Therefore, in the evaluation of the /3-function we expect that the deep infrared region 
is dominated by the smaller scale M2(0), so that the contribution of the larger scale Mi(0) 
can be neglected. Accordingly, we select //^ = M2{(j)) = JvX/(f) with t = ln(/i/A) and Eqs. 
( l3T]l32l) are replaced by 

| = -,„, + C.^«^ (42) 

^AMl_m . C^E^^iM^ . M!), (43, 

Consequently Eq. fl33l) becomes dr/dt = so that its solution in fl37j) now becomes r(t) = r, 
while the factor (A^ + 8) has to be replaced by (A^ - 1) in Eqs. fl3ll35ll36|) and Eq. §^ is 
formally unchanged. 

Since our analysis concerns the infrared region corresponding to large negative t, at this 
point it is essential to recall the different behavior of the theory in 2 < c? < 4 and d = 4, due 
to the presence, in the former case, of the Wilson-Fisher fixed point which merges with the 
Gaussian fixed point when d = 4. Therefore we have to treat the two cases separately and we 
start by considering 2 < d < 4. Then, from Eq. fH2l) one finds the Wilson-Fisher fixed point 



(that vanishes when e — )■ 0) and the scaling of u{t) and ^(t) near this fixed point is 

1 1/11 



1^^ 



u(t) u* \u u* 
15 



e'' (45) 



'1/2 



eW-(-j e-^\ (46) 

From the infrared behavior of the couphng in Eq. ( 145|) . one can write for the corresponding 
dimensionful quantity, after neglecting subleading terms: 

A 



\{t) 



u 

^ — e 
u* 



(47) 



Now we can focus on the inverse susceptibihty which, due to the rescahng of (p in Eq. ( I29l) . 
is XL^iJ) = v^ (dJ/dcp) = WI^Qj and it can be read from Eq. ( HTl) after inserting the exphcit 
form of the running parameters. Then, according to Eqs. ( H6l HTj), for small values of the 
source the linear term in J in the right hand side of Eq. ( HTl) can be neglected with respect 
to the other one, and one finds: 



x-l\J) 



XV' - — e 



eti2 




0A2 



(4-d)/2 



(48) 



In the limit J — > in Eq. (HHj) . the field </> can be replaced with 0o (up to subleading terms), so 
the inverse susceptibility vanishes as xl^i-J) °^ •^'^'^^ when J — )■ 0, which is the same behavior 
found in Eq. flT6|) in the limit of large A^. 

The comparison with the large A^ case can be pushed forward by looking at Eq. fj40|) and 
by including the scale dependence displayed in Eqs. (I44|45|46ll7r|) . For large N we get 



J 

7x 



6'" 



A 



'J\/a' 

0A2 



(4-d)/2 



U 

u 



6Ar+T 



3e 



' Jy/X 



{i-d)/2 



C,N 



(49) 



It is not difficult to see that Eq. (1491) coincides with the A^ — ?■ oo result in Eq. ( 1121) . at least 
in the region of small J and < e << 1, where one can neglect the linear term in p and 
approximate r(l — d/2) = r(— 1 + e/2) ^ — 2/e. Incidentally, we notice that the particular 
choice of the Goldstone mass as the infrared scale and of the /3- function in ( H2|) . has selected 
only the Goldstone bosons contribution, which is also dominant in the large A^ approximation. 
In this sense, it is not surprising to find the same result in Eqs. ( !T2|) and ( 149|) . 

Finally we analyze the problem in c? = 4, i.e. e = , where the coupling is dimensionless 
(A = u) and, as is well known and as observed in the N ^ oo case in Sec [TTl the divergences 
associated to the e~^ poles are typically replaced by logarithms of the Renormalization Group 
scale. In addition, the main difference with respect to the problem in (i < 4 is that the Wilson- 
Fisher and the Gaussian fixed points coincide in d = A, u* = 0, and this in turn implies that, 
when the ultraviolet cutoff is removed, one ends up with a noninteracting theory. In practice 
we can observe this 'triviality' property directly from the coupling constant, as obtained by 
integrating Eq. (j42j) in the limit e — )■ (below we define h 



N-l- 



167r2 



uit) 



u 



u 



1-ht h\t\ 



(50) 
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where the last term is obtained in the region of large negative t. In fact, from the second term 
in Eq. flSCT]) (i.e. before considering its approximation for large negative t that is displayed 
in the right hand side of Eq. ( !50l) ). it is evident that the limit of infinite ultraviolet cutoff 
A —7- oo, with fixed bare coupling u, corresponds to t — t- — oo for any finite infrared scale /x 
and therefore one finds u{t = — oo) = 0, i.e. the coupling vanishes at any finite fi. At the same 
time, we notice that the noninteracting limit u{t = — oo) = is recovered also at finite A, 
but for a vanishing infrared scale // = and, in conclusion, a nonvanishing coupling u{t) ^ 
requires finite values of t, i.e. finite infrared and ultraviolet scales. 

Then, going back to our problem, from Eqs. f H2|43|38p . in d = 4 and for large negative 
but finite t, we find r(t) = r and C,{t) ~ {h \t\) ' . With these results we can immediately 
speciahze Eq. (iO]) to rf = 4: 



\fu 6 



r(/i|t|)-' + ^ {h\t\) 
u 



(51) 



and, in the same way, the second derivative of the potential in Eq. fl4T|) becomes 

,-U)..K^„.^,-^. (52) 



where J has been inserted according to Eq. ( ISTj) . Eq. ( 1521) shows that the inverse longitudinal 
susceptibility, as obtained from Vrgi^ depends only on the ratio of the Renormalization Group 
scales A and /i through the logarithm t. Then, since the infrared scale is taken equal to the 
Goldstone boson mass, jj? = M2{(l)) = JvX/cj), we recover the logarithmic dependence of 
Xl i-J) o^ ^^^ source J (after neglecting the first term in the right hand side of Eq. (|52l) 
which vanishes linearly with J). It is also straightforward to realize that in the large N limit, 
Eq. (l52l) reproduces Eq. (ll7l) with the loop integral displayed in Eq. (!20|) . 

Therefore, our conclusion in d = 4 is that Xl^{J) vanishes as (ln(J))~^ when the external 
source is turned off, exactly as it happens in the large N limit. However it must be remarked 
that, because of the logarithmic dependence of our observable on the ratio of two scales, the 
limit J — )■ cannot be disentangled from the limit A — )■ oo and this means that the results 
obtained in the limit J — ?■ are those of a noninteracting theory. 

Before finishing we would like to investigate the possibility of observing a nonvanishing 
longitudinal susceptibility in d = 4. According to the above analysis, it is clear that formally 
the only way of achieving this result is to stay with a finite value of t, by keeping a finite 
A and also a finite infrared scale which for the moment could be taken as a generic value 
of the field: /i^ = 0^/m. If fi is taken much larger than Mi and M2 introduced above, then 
the problem related to the presence of two different infrared scales is no longer present and 
we can use the /3-functions in Eqs. (13T]) and (l32l) . This leads to a rescaling of the factor 
h-^ hr = h{N + 8)/{N - 1) in Eq. ([50]) and to the following explicit form of Eqs. ([37]) and 

r(t)~r(/i,|t|)^ ; ^(t) ~ (/^r 1^1)"'^' • (53) 
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As a consequence Eq. fllUl) now becomes 



\fu 6 



iV+2 



r(/i,|t|)-^ + ^ (/i,|t|)"' 



w 



(54) 



The minimum of Vrgi is located at 0o 7^ 0, such that J(0o) = in Eq. ([5 



2 ' '^^ 

1 2 



In 



iV+S 



(55) 



The square magnetization ^g/ti displayed in Eq. ([S5D, or equivalently the temperature r, 

corresponds to a finite quantity that can be used to set the infrared scale in this calculation. 
Then Eq. fHTl) . computed at = 0o (i.e. at J = 0) and with 00 replaced by r according to 
Eq. dSSD, gives (up to logarithmic corrections), 



Inf- 



rs^ 



Xl 2,2 



(56) 



Eq. ( 156|) provides a nonvanishing x^^, but it is evident that this result can only be obtained 
at the price of constraining t to a finite nonvanishing value, which in turn means that we 
are dealing with an interacting theory, effectively defined for values of the momentum smaller 
than A, and where the effects of the Goldstone bosons are screened by the finite infrared scale. 



V. CONCLUSIONS 

We studied the longitudinal susceptibility in the broken phase of the O(A^) theory, both in 
the simpler case of the large N limit and in the general case of finite A^. In this analysis we 
compared the results of three different approaches, namely the leading contribution of the 1/A^ 
expansion, the Functional Renormalization Group flow equations in the Local Potential ap- 
proximation and the improved effective potential via the Callan-Symanzik equations, properly 
extended to four dimensions through the expansion in powers of e. 

In the large A^ limit, the Functional Renormalization Group correctly reproduces a convex 
potential with a flat part in the region |$| < |<l>o|; while this does not occur in the potential 
obtained in the large A^ expansion. In fact, the latter potential is computed by performing an 
asjTiiptotic expansion of the path integral in Z\J\ around a single saddle point conflguration. 
In order to recover the convexity of the effective potential for |$| < |$o|; the contribution of 
the other saddle points has to be considered. The virtue of the Functional Renormalization 
Group is that it takes into account by construction all of these contributions. On the other 
hand, for |<I>| > |$o|) the agreement of the two approaches is excellent. 

According to the saddle point analysis performed in the large A^ limit , the inverse longi- 
tudinal susceptibility Xj} vanishes with the external source J as Xj} oc J*^/^, for 2 < c? < 4. 
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We have also shown that this scahng law is limited to a very narrow region around J = 0. In 
d = 4 the power law is replaced by the logarithm of the ratio of the ultraviolet cutoff and the 
infrared scale associated to the source J and xl^i-J) vanishes with J as (ln(J))~^. 

At finite A^, the Functional Renormalization Group numerical analysis shows the vanishing 
of Xl^ for d = 3 and d = 3.5. In d = A the region around the minimum of the potential is 
affected by strong numerical fluctuations and, as a consequence, it is not possible to conclude 
from the numerical analysis whether or not Xl^ has a discontinuity at |0o|- This point is 
clarified by means of the Callan-Symanzik approach. Within this latter framework we analyzed 
the longitudinal susceptibility for e > and recovered the same behavior already found in the 
limit N — )■ oo, namely Xl^ vanishing with J as J^^^ for e > and as (ln(J))~^ in li = 4, in 
agreement with the results of [g]. 

More generally, we find a full correspondence of our Renormalization Group improved 
expressions of the derivatives of the effective potential in the Callan-Symanzik approach when 
A^ — )■ oo, with the same quantities computed at the leading order in the large A^ expansion. 
Finally, in this approach it is clearly illustrated how the behavior of the inverse longitudinal 
susceptibility in d = 4 is simultaneously related to the vanishing of the Goldstone boson mass 
and to the limit A — )■ oo which leads to a noninteracting scalar theory. 
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